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Large optical gain from four-wave mixing instabilities in semicon- 
ductor quantum wells 
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PACS 71.35.-y Excitons and related phenomena 

PACS 78 . 47 . +p - Time-resolved optical spectroscopies and other ultrafast optical measurements 
in condensed matter 

PACS 42.65.Sf - Dynamics of nonlinear optical systems; optical instabilities, optical chaos and 
complexity, and optical spatio-temporal dynamics 

Abstract. - Based on a microscopic many-particle theory, we predict large optical gain in the 
probe and background-free four-wave mixing directions caused by excitonic instabilities in semi- 
conductor quantum wells. For a single quantum well with radiative-decay limited dephasing in 
a typical pump-probe setup we discuss the microscopic driving mechanisms and polarization and 
frequency dependence of these instabilities. 
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In gaseous atomic or molecular systems and simple 
Kerr media, four-wave mixing (FWM) processes lead, 
among other things, to transverse optical instabilities 
(e.g., refs. [1-5]). The interest in these FWM instabil- 
ities has recently been renewed by the demonstration of 
their effectiveness for all-optical switching at very low light 
intensities [6,7]. In this Letter we argue on the basis 
of a microscopic many-particle analysis that FWM in- 
stabilities can occur via nonlinear excitonic processes in 
a single semiconductor quantum well (QW). Instead of 
studying spontaneous off-axis pattern formation induced 
by these instabilities, we investigate their role in a pump- 
probe setup as illustrated in fig. 1. We find that the 
FWM instabilities can lead to large gain in the probe 
and (background-free) FWM directions that grows expo- 
nentially with the pump pulse duration, limited by the 
eventual buildup of incoherent exciton/biexciton densities. 
Our analysis shows that the materials conditions for ob- 
serving these instability-induced gains, though quite strin- 
gent, appear to be obtainable in currently available high- 
quality QW samples. 

FWM has been widely used in studies of microscopic 
processes in QWs (e.g., refs. [8-15]). An example of 
FWM driven instabilities in semiconductors is the para- 
metric amplification of exciton polaritons in planar QW 
microcavities [16-23]. In contrast to the microcavity, we 
focus on the "simple" system of a single QW. Using a 
microscopic many-particle theory we give a comprehen- 
sive stability analysis of the optical polarization field in- 
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Fig. 1: Illustration of the investigated pump-probe setup, in- 
cluding incoming and transmitted light pulses (reflected pulses 
are omitted for clarity) for probe, pump, and FWM direc- 
tion (the angle between the different propagation directions 
is strongly exaggerated). Selected contributions to the probe 
and FWM polarizations are schematically included in the fig- 
ure. The numbers in parentheses indicate the lowest order in 
the external fields in which they appear. 



duccd by a normally incident pump beam, treating all 
vectorial polarization state channels. In atomic systems, 
phase-space filling (PSF) nonlinearities drive the instabil- 
ities, and in microcavities it is the interaction between co- 
circularly (say, ++) polarized polaritons, related to the 
Hartree-Fock (HF) interaction and two-exciton (++) cor- 
relations [21-23]. We find that neither PSF nor excitonic 
interactions in the ++ channel are promising for instabili- 
ties in single QWs, mainly because of excessive excitation- 
induced dephasing (EID). Instead we show that instabili- 
ties and large probe gains can be produced in single QWs 
by virtual biexciton formation, which requires the pump to 
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contain both circular polarizations (e.g., to be linearly po- 
larized). Biexciton- induced bistabilities and instabilities 
have been studied before, but in bulk semiconductors (e.g., 
for CuCl [24-27]), rather than in quasi two-dimensional 
QWs studied here. These early works used bosonic mod- 
els and did not consider optical gain in the off-pump di- 
rections. It has been established (e.g., refs. [28-30]) that 
for typical QW systems (e.g., GaAs, ZnSe) a fermionic 
microscopic theory as used in this paper is needed for the 
understanding of excitonic PSF, spin-dependent exciton 
interactions giving rise to excitonic HF mean-field effects 
as well as retarded exciton interactions that can result in 
the bound two-exciton state (biexciton) and two-exciton 
continuum states, EID, and quantum memory effects. 

We investigate a pump-probe setup as illustrated in 
fig. 1 with the light propagating in quasi-normal incidence 
and with all optical pulses spectrally centered near the Is 
heavy-hole (hh) exciton resonance. Assuming all other 
resonances to be sufficiently far away, the coherent re- 
sponse of the system in the lowest-order nonlinear regime 
(x' 3 )-regime) has been well studied. We focus our analysis 
on small pump intensities where the many-particle effects 
listed above are dominant for the coherent optical QW 
response (i.e., we neglect higher than four-fermion or two- 
exciton correlations). The influence of higher-order corre- 
lations and incoherent exciton contributions is discussed 
later in our analysis. We start from the nonlinear equa- 
tion for the optically induced interband polarization p 
(+, — label the circular polarization states) and perform a 
spatial Fourier decomposition of p^ and the exciting field 
E^ with respect to the in-plane wave vector k. We la- 
bel the Fourier components with the subscripts s, / and 
p for probe (also called signal, with k = k s , assumed to 
be small but nonzero), background-free FWM (k = —k s ), 
and pump (k = 0) direction, respectively. The resulting 
equations are linearized in the weak probe field Ef- but 
solved self-consistently in the pump field E^ [31,32]. The 
equations for p J * read: 

see eq. (1) 

The nonlinear pump equation for p^ (not shown) involves 
the same nonlinear processes. We note that the solu- 
tion to eq. (1) goes beyond the x^-limit and includes 
the pump polarization up to arbitrary order. Here, e x 
is the ls-hh exciton energy, 7 a phenomenological exci- 
tonic dephasmg constant (not including the radiative de- 
cay), d cv the interband dipole matrix element, <^i s (r) the 
two-dimensional exciton wavefunction, v4 PSF the excitonic 
PSF constant, and V HF the excitonic HF matrix element. 
Unless otherwise noted, the time argument is t. The cor- 
relation kernels Q are given by Q ++ = Q = G + and 
G+- = g~+ = G+ +G~, with G ± as defined in eq. (22) 
in ref. [30], including a two-exciton dephasing rate 27. 
The propagation of the optical field E ± across the QW 
is described with a transfer-matrix method accounting for 
radiative corrections/decay and assuming the QW to be 
infinitely thin (e.g., ref. [33]). 
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Fig. 2: Transmitted intensity of a 500 fs probe (normalized to 
the peak intensity Iq of the incoming probe) . The shape of the 
20 ps pump pulse is included as the dotted line. Results are 
shown for XX pump-probe polarization, zero pump-probe de- 
lay, pump peak intensities 22kWcm -2 (solid) and 5.3kWcm -2 
(dashed-dotted), and equal pump and probe central frequen- 
cies tuned to — 0.6meV (solid) and — 1.2meV (dashed-dotted) 
below the exciton resonance, respectively. 

From a time-domain solution 1 of eq. (1) and the equa- 
tion for the nonlinear pump polarization, we find for co- 
linear (XX) pump-probe polarization configuration for cer- 
tain pump detunings Ae = hoj p — e x and above a certain 
pump threshold intensity / t h a- large energy transfer from 
the pump into the probe and FWM directions. To il- 
lustrate this phenomenon, fig. 2 shows the transmitted 
probe intensity for two different Ae which - after the ini- 
tial probe has passed - grows (almost) exponentially un- 
til the pump-induced QW polarization (not shown) falls 
below a certain threshold value (calculations for a high- 
quality GaAs QW at low temperatures where the radia- 
tive decay dominates the homogeneous linewidth as, e.g., 
in ref. [15]). For the results in fig. 2 about two orders 
of magnitude of time-integrated gain in the probe direc- 
tion is found. We find similar outgoing fluences in the 
background- free FWM directions (not shown). 

To analyze the gain mechanism, we supplement the 
full time-dependent solution by a linear stability analy- 
sis (LSA). The LSA is done without an incoming probe 
field and for a monochromatic cw pump field Ef(t) = 
E p ^e~ tu ' pt and pump polarization p^(t) = p^c~ luJpt , with 

Pp = Ep = (uip is the pump frequency). We evalu- 
ate the memory integrals in eq. (1) in Markov approx- 
imation [p s j(t') w p s ,f(t)e lu ' p ( t ~ t )] for the two-exciton 
continuum in the correlation kernels ^ ±:t , Q ^ . In con- 
trast, exact treatment of the bound biexciton state [in- 
cluded in the correlation kernel Q ±T in eq. (1) [30]] via 
the time-dependent amplitudes b s j(t) is crucial to include 
the related quantum memory effects (e.g., refs. [34-36]) 
which - as our analysis reveals - give rise to additional 
unstable modes. According to eq. (1) the term driving 
the bound biexciton amplitudes b s j(t) is proportional to 

1 Parameters: e x = 1.4965 eV, E£ x a 1.84 mcV, d cv = 4eA, 7 = 
0.01 meV, A PSF = 4agv / 2^/7, V HF = 2tt(1 - 3157r 2 /4096)ag 2 £*, 
with ag ss 170 A, Eg « 13meV. 
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PpPfj +PpPZ r With the ansatz p s j(t) = Psj(t)e 
and b s j(t) = b s j(t)e~ i2u3 ^ the probe and FWM dy- 
namics take the form ^p(^) = Mp(t), with p(t) = 

[p+(t),p+*(t),pj(t),pj*{t) 7 b s (t),b}(t)] T where the ma- 
trix M follows from eq. (1). If at least one of the eigen- 
values Xi of M fulfills Rc{A;} > 0, the system is unstable. 
An arbitrarily small seed of p s ^ would grow exponentially, 
until the matrix M ceases to describe the system correctly. 

On the PSF-lcvel [only line 1 of eq. (1)] the nonlin- 
ear terms from left to right are: (i) non-grating Pauli- 
blocking, (ii) self wave mixing (SWM) including power 
broadening with the general structure ip s j ~ PsjP p E p , 
and (iii) the phase-conjugate oscillation (PCO) with the 
structure ip s j ~ P*j s PpE p which gives the feedback be- 
tween probe and FWM signal (compare fig. 1) which drives 
the instability in atomic systems, e.g., in ref. [1]. Go- 
ing beyond the PSF-levcl, in the second line of eq. (1), 
HF exciton Coulomb interaction and excitonic correlations 
[the two-exciton continuum in the electron-spin triplet 
(-1 — h) channel] contribute to the two-exciton interaction. 
While the first two of these terms are of the SWM type 
(ipsj ~ PsjPpPp) and basically blue-shift the exciton res- 
onance in eq. (1) and lead to additional EID, the last two 
terms have the general structure of a PCO contribution 
(ipsj ~ P*f s PpPp) mediated by the two-exciton Coulomb 
interaction. The third line in eq. (1) contains the Coulomb 
interaction in the two-exciton electron-spin singlet (+— ) 
channel, including both the bound biexciton state and 
the two-exciton continuum in this channel, where the first 
term is of the SWM type and the last term has the struc- 
ture of a PCO coupling. 

We start the discussion with a a + polarized pump. 
Without Coulomb interaction, only the PSF-PCO can lead 
to an instability as it is the case in atomic systems. How- 
ever, within our model with zero in-plane exciton polari- 
ton dispersion, one can show analytically that PSF alone 
does not yield an instability (i.e., all Rc{Ai} < 0), which is 
analogous to the forward propagating case in atomic sys- 
tems [2]; PSF-SWM and PSF-PCO counteract each other. 
We note that an extension of Eq. (1) to realistic nearly 
quadratic dispersions [37] [e x = e x (k)] can - in princi- 
ple - yield instabilities for certain positive pump detuning 
£x(0) < fajj p < £ x (k a ), however our numerical studies have 
shown that they occur only at extremely high intensities, 
rendering the instabilities unlikely to be experimentally 
observable and justifying the neglect of the in-plane dis- 



persion. 

Without in-plane dispersion and including only HF 
terms (which, physically, leads to an exciton blue-shift 
caused by the HF-SWM terms), we found from a numer- 
ical diagonalization of M instabilities in a small spectral 
window for positive pump detuning, mediated by the HF- 
PCO contribution. In planar semiconductor microcavitics, 
the parametric polariton amplification has been attributed 
to this instability [16-22]. However, two-exciton correla- 
tions in the electron-spin triplet (++) channel [30] have 
been found to weaken this instability [23] , and in our case 
of a single QW we find that for realistic parameters the 
instability is inhibited when these correlations are taken 
into account; mainly because of EID due to two-exciton 
scattering. 

In contrast, a regime where excitonic correlations are 
strong but EID is kept small can be accessed for negative 
pump detuning Ae < with a linearly (say X) polar- 
ized pump that also excites the excitonic correlations in 
the electron-spin singlet (H — ) channel. The corresponding 
eigenvalues of M are shown in fig. 3 for a fixed total coher- 
ent exciton density defined as n^. otal = |p^| 2 + 2\b p \ 2 with 
the coherent exciton (|p^| 2 ) and bound biexciton (|6p| 2 ) 
densities (for Ae < we have checked that density con- 
tributions from the two-exciton continuum are negligible 
compared to n|j. otal ). The pump polarization p* induced 
by the X polarized pump follows from the (cubic) nonlin- 
ear pump equation and b p is the bound biexciton ampli- 
tude driven by a source term proportional to p p p p . 

In fig. 3(a) we find three different unstable regions 
(Re{A} > 0) caused by the biexcitonic (+— ) PCO terms 
in eq. (1). The labels XX and XY denote the vectorial 
polarizations of the pump (always X) and the unstable 
modes or probe fluctuations (either X or Y). The de- 
pendence of the eigenvalues with the largest real parts 
(which, when positive, are the maximum growth rates) on 
n totai j g s } lown m fig_ 5( a ) [ w ith fig. 5(b) relating n* otal 
to the pump intensity]. Clearly, the LSA helps us un- 
derstand the temporal growth in fig. 2, with maximum 
coherent densities n^ otal w 1.8 x 10 10 cm" 2 (22kWcm~ 2 ) 
and w 3.5 x 10 10 cm" 2 (5.3kWc m - 2 ). 

We consider a range of relatively small total coherent 
densities n!j. otal < 3 x 10 10 cm~ 2 to minimize the influ- 
ence of higher-order correlations not included in our the- 
ory. This restriction limits the maximum growth rates 
in the unstable regime [fig. 5(a)] and therefore dictates 
the minimum timescale that is required to achieve a de- 
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Fig. 4: (a) Pump intensity corresponding to the total coherent 
exciton density n^ otal = 1.6 x 10 10 cm~ 2 used in fig. 3. (b) 
Exciton (|pjf| 2 ) and biexciton (2|b p | 2 ) contributions to n^ otal . 



Fig. 3: Linear stability analysis for a linearly polarized pump 
for steady-state total coherent exciton density n^. otal = 1.6 x 
10 1() cm -2 . Shown are the real parts (a) and imaginary parts 
(b) of the eigenvalues Ai of the matrix M for negative pump 
detuning. Eigenvalues are represented by solid lines for the 
co-linear (XX) and by dashed-dotted lines for the cross-linear 
(XY) polarization configuration. The corresponding pump in- 
tensity and exciton (|pjf| 2 ) and biexciton (2|& p | 2 ) contributions 
to n* otal are depicted in fig. 4. The dotted line in panel (a) sep- 
arates the stable (Re{A} < 0) from the unstable (Re{A} > 0) 
regime. 

sired amount of gain. In principle, accumulation of an 
incoherent exciton/biexciton density in the system on this 
timcscalc (leading to additional EID) can cap the possi- 
ble amount of gain. However, for the time scales in fig. 2 
we believe that additional EID from incoherent densities 
would not crucially affect the presented results. 2 In addi- 
tion to the EID from incoherent densities, the exponential 
growth of gain with increasing pump length can also be 
limited by depletion of the pump signal, if, after a certain 
growth period, the probe and FWM signals have become 
sufficiently strong so that linearization of the polarization 
dynamics in these signals is no longer valid. Finally, in 
fig. 5(a) we show how the instability threshold density 
changes by increasing the effective dephasing in the sys- 
tem (e.g., by increasing temperature or due to EID). As 
a relatively sharp spectral feature around the two-photon 
biexciton resonance [fig. 3(a)] the XY instability may not 
be robust against even small deviations from the model 
conditions considered here. Therefore, in the following, 
we focus our attention on the unstable modes in the XX 
configuration. 

We can distinguish (i) an XX instability below (Ae < 

2 For pump excitation below the exciton resonance, energy conser- 
vation dictates that only scatterings from coherent and incoherent 
bound biexcitons (the incoherent biexciton density does not exceed 
2.5 X 10 9 cm -2 for the results shown in fig. 2) can efficiently cause 
EID losses for the coherent exciton and biexciton densities. As- 
suming that the scattering rates for processes involving at least one 
bound biexciton are considerably lower than for two-exciton scat- 
tering, these EID contributions are negligible for the time scales 
considered in this paper. 



-E% x /2) and (ii) above (Ae > -E xx /2) the two-photon 
biexciton resonance at Ae = — Ef x /2 s» — 0.92meV. Fo- 
cusing on (i), we see from fig. 5(b) that it occurs in a 
region of bistability (i.e., three possible steady state co- 
herent exciton densities for fixed intensity). In the cor- 
responding time-domain calculation in fig. 2 the instabil- 
ity of the intermediate branch of the bistable pump solu- 
tion (cf. ref. [27]) is reflected in the exponentially growing 
probe signal over time, and as demonstrated in fig. 2 can 
give rise to large probe gain. However, we note that in this 
bistable regime the system dynamics sensitively depends 
on the 'initial conditions', determined by the overall tem- 
poral pump pulse shape. In fig. 2 we use a Gaussian pump 
pulse with intensity and pulse length chosen such that the 
pulse induces a coherent density on the intermediate (un- 
stable) branch in the bistable regime [cf. dashed-dotted 
line in fig. 5(b)]. However, for a different choice of pump 
parameters the induced polarization may have mainly con- 
tributions from the lower and upper (stable) branches, so 
that no exponential growth (instability) in fig. 2 would 
be observed. Turning to the instability (ii), we find a 
much lower threshold (as low as n^. otal w 3 x 10 9 cm~ 2 ) 
than for the other configurations. Furthermore, fig. 3(b) 
explains the apparent beating on the probe signal (solid 
line in fig. 2) during the exponential growth period. In 
case (ii), the unstable modes exist in pairs with complex 
conjugate eigenvalues. Whereas the growth rates, Re{A}, 
are the same for the two modes in each pair, the frequen- 
cies, Im{A}, have opposite signs. The unstable growth of 
the probe and FWM signals occurs at the two frequen- 
cies ujp ± |Im{A}|. This is inherently related to a non- 
Markovian quantum memory effect, brought about by the 
time-retarded structure of the biexcitonic SWM term in 
eq- (1). 

Our analysis is applicable to other semiconductor ma- 
terials provided suitable scalings in lengths (with ag) and 
energies (with E^) are made. Since we believe that FWM 
instabilities in QWs can be expected as long as the exci- 
tonic dephasing is not much larger than the radiative de- 
cay (e.g., w 0.05 meV for GaAs [15], « 0.3 meV for ZnSe), 
materials with a very small ratio j/E% would be ideal. 

In conclusion, we predict optical FWM instabilities and 
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Fig. 5: (a) Real part of the eigenvalue Ai of M with the largest 
real part corresponding to co- linear (XX) and cross-linear (XY) 
pump-probe polarization configuration, respectively, (b) Total 
coherent exciton density n^° tal vs. pump intensity inside the 
QW. (a), (b) Results are shown for different pump detunings, 
Ae = -0.6 meV (solid), Ae = -0.92 meV (dashed), and Ae = 
— f .2 meV (dashed-dotted). The thin solid line in (a) shows the 
same result as the solid line but for a ten times larger dephasing 
7 = 0.1 meV. 



corresponding strong probe and FWM signal gain in a 
semiconductor QW. Further research, e.g., into possible 
transverse pattern formation and propagation effects in 
multiple QWs, would be desirable. 
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